Spin dynamics and structure formation in a spin-1 condensate in a magnetic field 
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We study the dynamics of a trapped spin-1 condensate in a magnetic field. First, we analyze 
the homogeneous system, for which the dynamics can be understood in terms of orbits in phase 
space. We analytically solve for the dynamical evolution of the populations of the various Zeeman 
components of the homogeneous system. This result is then applied via a local density approximation 
to trapped quasi-lD condensates. Our analysis of the trapped system in a magnetic field shows that 
both the mean-field and Zeeman regimes are simultaneously realized, and we argue that the border 
between these two regions is where spin domains and phase defects are generated. We propose a 
method to experimentally tune the position of this border. 
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I, INTRODUCTION 

Bose-Einstein condensates (BECs) with a spin degree 
of freedom are an interesting field of research in many- 
body physics as they realize both superfluidity and mag- 
netism in a well-controlled environment. First realized 
experimentally with 23 Na ten years ago their study 

has matured remarkably over the last few years, with sev- 
eral groups studying their dynamics 0, 0, H, @] and ther- 
modynamics @, S|. Of particular interest is the study 
of the process by which spin domains are formed dur- 
ing time evolution, a phenomenon observed experimen- 
tally @, S [Io| and in numerical simulations based on a 
mean-field approach (ill . H3 . [T3 | . 

The complicated dynamics of these non-linear systems, 
especially when they are subjected to time- varying exter- 
nal fields, makes the physical understanding of the struc- 
ture formation process somehow elusive. To address this 
point, we present here a simple model based on an ana- 
lytic solution for the homogeneous system for arbitrary 
magnetic fields B and magnetizations M. This solution 
is then applied to the study of realistic, trapped spin-1 
condensates by means of the local density approximation 
(LDA) . This approximation has already been applied suc- 
cessfully in a number of studies on scalar BECs, as well as 
cold Fermi gases. From the analysis of our results we are 
able to provide an intuitive picture of the process lead- 
ing to the structure formation. Further, we argue that 
it should be possible to experimentally "tune" the spatial 
region where this process starts within the condensate. 

The paper is organized as follows. In Sect. Ill Al we 
present the phase space of a homogeneous system under a 
magnetic field B and for arbitrary M , and introduce the 
phase-space orbits that describe the dynamics of a con- 
servative system. In Sect. Ill Bl we solve analytically the 
dynamical evolution of the homogeneous system. Then, 
in Sect.Hll]we describe our local-density approximation 
for a trapped system and present numerical results for 



its dynamics (Sect. llTTA"]) . which we compare with simu- 
lations based on a mean-field treatment (Sect. UlTCj) . In 
Sect. IIVI we discuss the progressive dephasing of differ- 
ent spatial points of the condensate in a homogeneous 
magnetic field, and relate this to the process of structure 
formation, with an indication of a possible experimental 
test. Finally, we conclude in Sect. fVl 



II. ANALYTICAL RESULTS FOR THE 
HOMOGENEOUS SYSTEM 

A. Energetics of the homogeneous system 

A homogeneous condensate of atoms with total spin 
F can be described by a vector order parameter tp with 
2F + 1 components, 



1> 



(i) 



The density of atoms in a given Zeeman component m = 
—F,--- ,F is n m = | ip m | 2 and the total density is given 
by n — J2 m IVVnl 2 - Introducing the relative densities for 
the homogeneous system A m = n m /n, one has 



i 



(2) 



Given that n is a conserved quantity, Eq. (J2j) will be ful- 
filled at all times during the dynamical evolution. More- 
over, the magnetization 



M = mA„ 



(3) 
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is also a conserved quantity [1 11 ] - 

We now focus our analysis to the case of a F — 1 con- 
densate. We write the various components of the order 
parameter as ip m — \Jn\ m exp(z0 m ). This ansatz, to- 
gether with conditions |(2|) and |(3]) , leads to the following 
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M=0.3, B=100 mG, n=4.4x10 14 cm" 3 
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Figure 1: Energy (in units of |c2|n) of the homogeneous sys- 
tem for the cases (a) M = 0, B = and (b) M = 0.3, B = 
1 mG, as given by Eq. (j4j for a spin-1 condensate of 87 Rb. 



expression for the energy per particle of the homogeneous 
system in the mean-field approach [3, [TB|: 



M 2 



-371/2 



Figure 2: Contour plot of the energy surface corresponding 
to M — 0.3 and B = 100 mG. The white line shows the orbit 
corresponding to the initial conditions AfJ 1 = 1/2, 6i n = 7r/2 
(indicated by the white dot). The minimum of £ is at Ao ~ 
0.455, 6 = 0. Note the presence of open orbits for energies 
above that of the indicated white line. 



B. Dynamics of the homogeneous system 

We are interested in the time evolution of the densities 
of the different Zeeman components, n m = n\ m . From 
Eqs. ([2l) and J3| we have that 



A±i - 



1 ± M - A 



(5) 



£(X ,M, 0) = c 2 n A (l - A ) + — + 



Therefore, we only need to follow the evolution of Ao, 
which is given by 



= Ao = ^c 2 n\ a ^/(l - A ) 2 - M 2 sm0 . 



+ A o V(l-A o ) 2 -M 2 cos0 + 5(1-A ). (4) 

Here 6 = 28q — 8± — 6-\, while c 2 is given in terms of 
the s-wave scattering lengths a/ in the channels of to- 
tal spin / = 0,2, by c 2 = 47rfr 2 (a 2 - a )/(3M), with M With Eq flU, we rewrite this as 
as the atomic mass. Finally, 5 — {E- + E + — 2E a )/2, 
where the energies of the atomic Zeeman states are 
given by the Breit-Rabi formula E m — — i?hts/8 — 
Eh&y/l + ma + a 2 /2 (m = — 1,0, +1), with E h(s be- 
ing the atomic hyperfine splitting and a = {gi^N + 
9j^b)B I 'Ehf s is a function of the external magnetic field 
B. Here, gi, gj are the nuclear and electronic Lande fac- 
tors, and [in,Ii<b are the nuclear and Bohr magnetons, 
respectively. A sketch of the surface £ is given in Fig. [TJ 



(6) 



( A °) 2 = -^\(^n\ ) 2 [(1-Ao) 2 -^ 2 



B - 5(1 - Ao) - can A (l-A ) + 



It can be shown that the term in Aq actually drops out 
and we are left with a cubic polynomial on Ao, 



(Ao) 2 = A(X - Ar)(A - A 2 )(A - A 3 ) , 



As indicated above, M. is a constant during dynamical 
evolution. Similarly, given initial conditions (A n ,#i n ), 
E = £(A n n , M, #i n ) will also be conserved, thus defining 
an orbit on the surface £ in (Ao, 0) space. A sketch of one 
such orbit is presented in Fig. [2j One should note that, 
depending on the initial conditions, the orbit defined by 
E = const can be closed or open. In the first case, 6 — 
6(t) will be a periodic function of time, while in the latter 
case, \0(t)\ will grow indefinitely with time. In both cases, 
however, Aq = Ao(t) will be a periodic function of time. 



with 



A:=- 



Bc 2 n5 



(7) 



(8) 



and Aj (j = 1, 2, 3) are the roots of (A ) 2 , Ai < A 2 < A 3 . 
For ground state (F = 1) alkalies 5 > 0. Therefore, 
Ai < AqW < A 2 for c 2 < and A 2 < X (t) < A 3 for 
c 2 > For concreteness, in the following we will 

assume c 2 < 0, i.e., ferromagnetic interactions. 
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We will now integrate the time evolution of Ao. To 
do so, we introduce an auxiliary variable z through 
Ao = (A 2 - Ai)z 2 + Ai. This will satisfy the differen- 
tial equation 



v/A 2 -A!vV-l)(z 2 -fc- 2 ) 



where we defined 



(9) 



(10) 



The first order differential equation |(9|) can be solved 
analytically by separating the variables z and t, and in- 
tegrating: 



dz 



A r i 

2 J to \/A 2 - Ai j. 



1 



V(l-^ 2 )(^ 2 
dz 



z 2 ) 



VA 3 - Ai J zo V(l-^ 2 )(1-^ 2 ) 

The solution to the last integral can be expressed in terms 
of the elliptic integral of the first kind [24j , 



F(0,fc) 



dz 



v /(l-z 2 )(l-fc 2 z 2 ) 



Taking as initial condition z(t = to) = zq and using 
the fact that F(—u,k) = —F(u,k), we can express z t 
in a compact form by means of the Jacobi elliptic func- 
tions which are defined as the inverses of the elliptic 
integrals, 



z t = z(t) 



sn 



( 70 + ^pS (t - t „)| t ) <„, 



with zq = sn(7o|fe), i.e., 70 := F(arcsin(zo), k). Finally, 
we undo the change in variables to write down the time 
evolution of the population of the \m — 0) state, 



Ao(t) - Ai+ 



(A 2 -A 1 ) B n 2 U+ ^- Al) ft-t ) 



In accordance with the identity [18j 

1 - cn(2a|fc) 



(12) 



sn (a\k) = 



l + dn(2a|fc) ' 



and given that both cn(2a|fe) and dn(2a|fe) are periodic 
functions in a with period 2K(k), Xo(t) will be a periodic 
function of time with period 



T = 



2h 



x /-2c 2 n<5(A 3 - Ai) 



K 



/A2-A1 
A 3 - A x 



(13) 



Here, K(k) = F(ir/2,k) stands for the complete elliptic 
integral of the first kind. We note that result (fl~3|) agrees 



with that in Ref. [17|, where T was calculated directly 
by performing the integral T = <j> dAo/Ao over a period of 
evolution. Further, let us point out that the average value 
Aq v = (1/T) J t * ° +T Ao(£) dt does not necessarily coincide 
with the position of the minimum of £, i.e., Ag v may differ 
from the equilibrium value Ap q (as given, e.g., in Ref. 
for the case B — 0). This is illustrated in Figs. [3]Jb) 
andE3Jb). 



C. Evolution in the absence of a magnetic field 

We observe that the representation of (Ao) 2 as a cubic 
polynomial on Ao, Eq. cannot be performed when 
A = 0, i.e., when B = 0. In this case, the analytic expres- 
sion (fl2| is meaningless, as it would apparently result in 
no time evolution at all. Actually, in this situation, (Ao) 2 
can be written as a quadratic polynomial on Aq: 



(Ao) 2 = {2c 2 nEX 2 - 2c 2 nE'X + {E'f} 
= Z(Ao-A 1 )(A -A 2 ) 



Al ' 2 = 2E 



2E 

c 2 n 



(Ai < A 2 



(14) 
(15) 



Here E' = E - c 2 nM 2 /2. Note that A := -8c 2 nE/h 2 oc 
— c| < for c 2 < as well as for c 2 > 0, and in both cases 
we will have Ai < \o(t) < A 2 . Following a procedure 
analogous to that above, we arrive at 



A (i)=A 1 + (A 2 -Ai)sin 2 7 



Tb=o 



(t-to) (16) 



with sin7 = [(A n - Ai)/(A 2 - Ai)] . In this case, A 
follows a pure sinusoidal evolution as has been predicted 
before in a number of references, e.g., [3, E3, El- The 
average value is Ag v = (Ai + A 2 )/2 = E'/(2E), and the 
period reads (compare with [19]) 



yJ2c 2 nE 



(17) 



We show in Fig. [3] the time evolution of Ao(i) for two 
representative cases. The different panels compare the 
analytic evolution -given by Eq. I112II or lfT6|) - with a 
numerical solution of the corresponding equation for Ao- 
In all cases, we see that the amplitude as well as the 
period of the time evolution are well predicted by the 
analytic results. Finally, we show in Fig. |4]a plot of 9{t) 
vs. Ao(t) corresponding to the time evolution depicted in 
Fig. [3l For the case with magnetic field and M. ^ we 
observe that the average value Ag v w 0.433 (indicated by 
the dashed line) differs from the position of the minimum 
of £ (Ap q w 0.455, cf. Fig. [2]) due to the deformation of 
the orbit. 
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Figure 3: Evolution of the population of the m = Zeeman 
component, Ao(t), for the cases (a) M = 0, B = 0, and (b) 
M — 0.3, B — 100 mG, starting in both instances from A n = 
0.5, 6i n = 7r/2. In both panels, the solid line corresponds to 
the analytic result, (a) Eq. (1 161) or (b) (11211 . while the circles 
are a numerical integration of the differential equation for 
Ao. The dashed line gives the expected average value of Ao, 
Ag v » 0.433. The arrows indicate the amplitude and period 
as predicted by the analytical results. In the bottom plot, 
also the value of the equilibrium population Ag q « 0.455 is 
indicated by a dotted line, while the dashed-dotted line stands 
for a fit to Eq. l121[| (displaced vertically by 0.1 for clarity). 



<x> 




CD 



Figure 4: Phase-space plot (Xq,6) corresponding to (a) the 
evolution shown in Fig.^a) and (b) Fig.[3][b) (compare with 
Fig. [2j. The solid line is the analytic result in (a) Eq. lllCIl 
and (b) Eq. lfl2|l . while the circles are the solution of the 
differential equations for Ao and 9. The vertical dashed lines 
stand for the average value Ag v in each case. 



III. DYNAMICS OF THE TRAPPED SYSTEM 



We have established in the previous section the dy- 
namical evolution of a homogeneous spin-1 condensate, in 
terms of orbits in the (Ao, 0) plane constrained by (i) con- 
servation of density, (ii), conservation of magnetization, 
and (iii) conservation of energy. The resulting dynam- 
ics of the population of the \m = 0) Zeeman component 
has been shown to be a periodic function of time, with 
a period determined by the density n of the system, its 
magnetization Ai, as well as the initial conditions of the 
evolution (implicit in E and, therefore, in {Aj}j = i i 2,3 or 
{Aj}j = i, 2 ), cf. Eqs. HU) and (fT7|). Now, we will transfer 
these results to a realistic case of a trapped, quasi-one- 
dimensional (ID) condensate. 



A. Local-density approximation 

The initial conditions for the evolution of a trapped 
spinor condensate are the set of complex values ip™( x ) 
for all Zeeman components m and all positions x where 
the density is not zero. In typical experiments, the prepa- 
ration of the initial state is such that n m (x)/n(x) is a 
constant independent of position. This, together with 
the fact that |ca | <C |co| for the systems studied so far, 
has lead to some theoretical works based on the so-called 
single-mode approximation (SMA), which assumes that 
n m (x,t) — n(x)X m (t) for all times t of the evolution, 
i.e., that the spatial variation in the density of each 
Zeeman component is always given by the total density 
profile. However, numerical studies beyond the SMA 
(e.g., [ll|, [12, Gil) predicted the formation of spin do- 
mains as time goes by. These have been observed in a 
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number of experiments, e.g., @,@|- In order to be able to 
observe the formation of spin domains during time evolu- 
tion in a trapped system, we will therefore not make use 
of the SMA, but apply the analytical results of Sect. |H] 
via the LDA, i.e., we will assume that the evolution of 
the \m = 0) population at each point within the conden- 
sate, Xo(x,t), is given by Eq. (fl2|) [or Eq. lfl6|) ] with the 
substitution n — > n(x). Here, the total density is nor- 
malized to the total number of atoms in the condensate, 
J dxn(x) = N. Similarly, we introduce the local densi- 
ties of atoms in a given Zeeman state n m (x, t) normalized 
as J dxn m (x,t) = N m (t). The conservation laws read 
now J2 m N m(t) = N and E„mJV m (i) = M. We note 
that n(x) does not change in time at low enough temper- 
atures [llfl unless momentum is imparted to the center of 
mass or to one or more of the Zeeman components [2p| . 

In the language of the phase space introduced in 
Sect. Ill A} a trapped system corresponds to an infinite- 
dimensional phase space, with a pair of variables 
(Xo(x),6(x)) associated to each point x. According to 
the LDA, we divide this whole phase space in sections 
corresponding to the different positions, and assume that 
they are independent. The initial condition described 
above, n m (x)/n(x) = const, corresponds then to the 
dynamical system starting in all the different positions 
x at the same point of the corresponding phase space, 
(Xo(x,t = 0) = A n , 8(x 7 t = 0) = 6 in ). The dynamical 
evolution of the system corresponds then to the popula- 
tion Xo(x 7 t) at each point x following its own particular 
orbit in the corresponding (Xo(x),6(x)) space, that is, 
Ao at position x follows the dynamical equation of the 
homogeneous system (fl2"l) [or Eq. (fl6j) ] with the param- 
eters Aj and A determined by the local density n(x). In 
other words, we assume that the position dependence is 
only parametric, and comes through the values of the pa- 
rameters Aj = Aj(x) and T = T{x). We will indicate 

this by Xo(x,t) = X^ BA {x,t) = A„ (x) (f). The density at 
position x of atoms in the Zeeman component m at time 
f will then be 



approximation, 



n m (x,t) = n(x)X m (x,t) 



(18) 



with X±i(x,t) = X±i\t) given by Eq. J5j with the sub- 
stitution Ao — > X^ x \t), and M = M{t = 0) is a con- 
served quantity [111 ] - 

Note that the orbits associated to different points x 
may differ from one another, as their shapes depend in- 
ter alia on the local density n(x), cf. Eq. Q. This fact, 
together with the position dependence of the parameters 
Aj(x) and T(x), is expected to lead to a dephasing of the 
evolution of the partial densities n m (x,t) at the different 
points, washing out the oscillations in the integrated pop- 
ulations, N m (t), in contrast to the stable oscillations that 
we have found for the homogeneous system, cf. Fig. [3j 

In order to evaluate N m (t) it is necessary to know the 
density profile of the system. A good estimate for n(x) 
in trapped atomic gases is given by the Thomas-Fermi 



n TF (x) 



- ) R\ 
0, 



-(R^-lxl 2 ), 



\x\ < R TF 
otherwise 



(19) 



For a quasi-lD system with total number of atoms N 
and central density n max , Rtf = 3N / (4n max ) . The inte- 
grated population in \m = 0) then reads 



N (t) = 



dxn TF (x)Xo TF[x) (t) 



(20) 



B. Analytic approximation with sinusoidal time 
dependence 

The time dependence of Aq (t) has in principle to be 
calculated from Eq. fT2|) for each position x at each time 
step, and then the integral (|20| performed numerically to 
determine N (f ) . It is possible however to give an analyt- 
ical estimation for N (f ) if we make a further assumption 
on the time evolution. From Fig. [3j we see that the evo- 
lution of A (t) for the homogeneous system is very close 
to a sinusoidal function even when B ^ This is 

illustrated in Fig. [3jb) , where a function of the form 



X c os {t) =o + 6cos(7 + i4) 



(21) 



has been fitted to the numerical values obtained from 
Eq. l[T2|) . The fit is very good, even for this case, 
where the orbit in phase space is strongly deformed 
[cf. Fig. Hjb)]. The advantage of approximating the time 
evolution of Ao by Eq. (|2~T]l is that it allows for an analytic 
evaluation of the spatial integral ([20)l , taking into account 
the position dependence of v. Indeed, from Eq. ([T3|l we 
expect v(x) cx n(x) cx (R^ F — x 2 ). It is easy to show that 



N™ s (t) = / dxn TF (x) 



a + b cos 



7 + ^o(l - sr-)* 



R 2 



la(v t) 3/2 + 6b^/l^tsm(2'y) 



6M) 3 / 2 

+ 3\/2^ : 6{cos(^)S'(ry) - sin(^)C(ry)} 

+ 6V2^bv t {cos(^)C*(?7) + sin (n) S (rj)} . (22) 

Here S (r/) and C (77) are the Fresnel integrals [l8| , and we 
introduced fi = 7 + v§t and 77 = ^/2i/ot/n. 

We show in Fig. [5] the time evolution of the integrated 
|m = 0) population as given by Eqs. f2"0|) and (|22| . 
This calculation has been done for a quasi-lD system 
of 20 000 87 Rb atoms in a trap such that the central 
density is 4.4 x 10 14 cm -3 . The initial conditions are 
A n = 0.5, 6>j n = 7r/2 and M. — 0.3 and we have taken a 
magnetic field B — 100 mG (cf. Fig. [2$. The solid line 
in the figure corresponds to the numerical integration 
of Eq. (22) with A (x,i) = X^ (x] '(f) given by Eq. (fl2f . 
The circles stand for the analytic expression (|22|) with 
the parameters a, 6, 7,^0 taken so that Ao° s (f) for a ho- 
mogeneous system with density n = n max reproduces the 
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Figure 5: Time evolution of the integrated |m = 0) popu- 
lation (normalized to the total population, N = 20 000) for 
Ajj 1 = 0.5, 9 iB = tt/2, M = 0.3 and B = 100 mG. The solid 



line shows the LDA result, Eq. (20J), with Xo(x,t) = Xq {x] '(*) 



The circles stand for the analytic estimate of Eq. 11221 . and 
the dashed line is the result of integrating the set of coupled 
Gross-Pitaevskii equations lf24|) . 



same behavior as that given by Eq. (fl2|) at the same den- 
sity: a = (Ai + A 2 )/2, b = (Ai - A 2 )/2, 7 = 2 7o , 14, = 
27t/Tlda, and Tlda = T(n max ). The agreement between 
the two calculations is very good at all times. There- 
fore, we conclude that the average value of Ao as well as 
the characteristic period of the oscillations is well deter- 
mined by the values Aj and Tlda calculated with the 
central density, while the time scale for the damping of 
the oscillations is determined by the spatial profile of the 
density. 

Regarding the dephasing of the evolution of Xo(x,t) 
among different points, it is not very strong, in the sense 
that the damping of the oscillations is relatively slow. To 
be more precise, one can have a reasonable fit to the solid 
line in Fig. [5] by a function of the form 

N (t) = N + AN exp(-aVt) cos f 2 7o + %t J (23) 



with No 



ANn 



V 



b , 7o ~ 7o and T w T LDA . 



C. Comparison with the mean-field approach 

We proceed finally to compare the approximate calcu- 
lation of No(t) with a more complete approach in terms 
of the dynamical equations for the three components of 
the vector order parameter, ip m (x,t), cf. Eq. (U). In the 
mean-field approximation, such equations can be cast in 
the form of three coupled Gross-Pitaevskii equations [ll|, 

(24a) 

ih^- =[Hs + c 2 (ni + + 202^1^0^-1 , (24b) 



where H s = -h 2 /(2M) V 2 + V ext (x) + c Q n(x) and c = 
47rfi 2 (a + 2a 2 )/(3M). 

The results of solving Eqs. lf24j) with a Runge-Kutta 
algorithm are included in Fig. [5] as a dashed line. The 
average value of the oscillating Xo(t) is well estimated by 
the analytical model of Sect . Ill Bl Also, the characteristic 
time scale of the oscillations is well estimated by Eq. (fl3|) . 
The overall agreement is good for times t < 100 ms. 
After this time, the analytical estimate keeps oscillating 
with a slowly decreasing amplitude, while the numerical 
solution of the coupled equations ([24| shows fluctuating 
oscillations. This behavior has been observed before, and 
the transition at t — t- mst ~ 100 ms has been related to 
a dynamical instability that leads to the formation of 
dynamical spin domains in the system 0, [ill, G3 • It is 
thus not surprising that our simple model fails for t > 
f inat • It is nevertheless remarkable that the time scale set 
by ?lda = T(n max ) s» 89 ms is still a good estimate of 
the characteristic oscillation time even much later during 
the time evolution. 



IV. DEPHASING IN A MAGNETIC FIELD AND 
THE PROCESS OF STRUCTURE FORMATION 
IN FINITE SYSTEMS 

A qualitative difference between the homogeneous sys- 
tem and the confined one appears when a magnetic field 
is present and, therefore, A ^ 0. The dynamics of a 
spinor condensate in a magnetic field is known to show 
two limiting behaviors: the mean-field regime, where the 
interaction energy dominates the evolution, and the Zee- 
man regime, where the evolution is driven by the Zee- 
man term of the Hamiltonian E3, HH • The crossover 
between the two regimes occurs when c 2 ?i ~ S. This 
transition can be studied in real time by changing the 
(homogeneous) mag netic field on which the condensate 
is immersed JlO, Ea E3| ■ 

This transition can also be observed between different 
spatial regions of an inhomogeneous system. Indeed, if 
we assume that the magnetic field, magnetization and 
central density are chosen so that |c 2 |n(a; = 0) > S (so 
that at the center we are in the mean-field regime), then 
at the wings of the system, where n{x) — > 0, we will be 
in the Zeeman regime. Therefore, we expect to have a 
region in real space where the behavior with time changes 
qualitatively. For a profile as in Eq. (fT9| . this transition 
border is given by 



R 



= Re 



TF 



c 2 n n 



1/2 



(25) 



Naturally, for S = 0, there is no transition (the density 
vanishes at x = Rtf)- On the other hand, for large 
enough magnetic field the whole system is in the Zeeman 
regime (^trans = 0). 

These two regimes evolve with different characteris- 
tic times, T m{ ~ H/(\c 2 \n) and T Zeeman ~ H/y/2c 2 nS, 
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cf. Eqs. (fT7|) and (fT3|) , Because of this, we can expect 
Ao and the phase in the inner part of the condensate 
< £ trans ) to evolve at a different rate than in the 
outer wings of the system (|x| > a; tr ans), resulting in 
a particular spatial dependence of the phase. We note 
that the appearance of a spatial structure in the phase 
will lead to the creation of spin currents [ll[ and, thus, 
to spin textures as reported in [3, [13]. Even though a 
smooth density profile will lead to a smooth variation in 
T(x) — T(n(x)) with position, from our model we ex- 
pect that these qualitatively different behaviors should 
be observable for times t > t tians = min{T mf , T Zee man } 

Interestingly, in light of the discussion in Sect. IIII C| 
we observe that the time when the dynamical instability 
is expected to set in is close to the time when the di- 
vergence between mean-field and Zeeman regimes should 
be observable: t; nst ~ itrans- Because processes such as 
spin currents fall beyond LDA, their appearance implies 
a breakdown of our model, which is therefore not appli- 
cable to analyze the process of structure formation. This 
breakdown explains the lack of agreement between the 
results of our LDA model and those from Eqs. (|24j) for 
t ^ iinst observed in Fig. [5] 

The experiments reported in Ref. [13] showed the 
appearance of spin domains to be simultaneous with 
that of topological defects (phase windings) and also 
spin currents. This observation is consistent with the 
model just sketched. The time scale for the appear- 
ance of spin domains is estimated in that reference to 
be ~ h/(2\c 2 \n) [13]. Similarly, Saito et at. Q deter- 
mined the time scale for the occurrence of a dynamical 
instability to be t- lnst = hj (|c2|n) when the magnetic field 
is small; this estimate coincides with our T m f. On the 
other hand, for larger magnetic fields [q > \c2\n with 
q = {jibB) 2 I '(4£ ; i 1 fs)], the relevant instability time scale 
is i inst = h/y/q\q + 2c 2 n\, which is similar to Tzeeman- 

From their simulations, Saito and Ueda indicated [13] 
that the formation of spin domains starts at the center 
of the condensate, and then spreads out. In our model, 
however, the position where the phase slip appears is de- 
termined by Xt ra ns ; and therefore is in principle amenable 



to be modified experimentally. It seems interesting to in- 
vestigate the prospect to control the spatial appearance 
of spin domains and phase structures as predicted by 
Eq. HMD- 



V. SUMMARY AND CONCLUSIONS 

We have studied the dynamics of a trapped spin-1 con- 
densate under a magnetic field. First, we have analyzed 
the homogeneous system and seen that its dynamics can 
be understood in terms of orbits in the (Ao, 6) space. We 
have then solved analytically for the dynamical evolu- 
tion Ao(t). We have used this information to study the 
trapped system by means of the Local Density Approx- 
imation (LDA). The results of this approach agree with 
those of the mean-field treatment for evolution times be- 
fore the occurrence of a dynamical instability [13]. In 
particular, the expected average value of Ao, as well as 
the characteristic time scale of its dynamics are well pre- 
dicted by the formulas for the homogeneous system. 

Our analysis of the trapped system has shown that, in 
the presence of a magnetic field, both the mean-field and 
Zeeman regimes are realized in a single spinor conden- 
sate. The analysis of this model allows for some qual- 
itative insight into the process of structure formation. 
In particular, our model identifies a transition point [cf. 
Eq. |25]) ] around which this structure is generated, and 
predicts that it should be tunable, which could be tested 
in future experiments. 
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